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Since the standard inflationary paradigm is based on quantum field theory on classical space- 
times, it excludes the Planck era. Using techniques from loop quantum gravity, the paradigm is 
extended to a self-consistent theory from the Planck scale to the onset of slow roll inflation, covering 
some 11 orders of magnitude in energy density and curvature. This pre- inflationary dynamics also 
opens a small window for novel effects, e.g. a source for non-Gaussianities, which could extend the 
reach of cosmological observations to the deep Planck regime of the early universe. 
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The inflationary paradigm has had remarkable success 
in accounting for the inhomogeneities in the cosmic mi- 
crowave background (CMB) that serve as seeds for the 
large scale structure of the universe. However it has 
certain conceptual limitations from particle physics as 
well as quantum gravity perspectives. For example: i) 
The physical origin of the inflaton and its properties re- 
mains unclear; ii) Since the background geometry and 
matter satisfy Einstein's equations, the big bang singu- 
larity persists [1]; iii) One ignores pre- inflationary dy- 
namics and simply requires that perturbations be in the 
Bunch Davies (BD) vacuum at the onset of the slow roll; 
and, iv) When evolved back in time these perturbative 
modes acquire trans-Planckian frequencies and the un- 
derlying framework of quantum field theory on classical 
space-times becomes unreliable. Here we will not address 
any of the particle physics issues. Rather, we focus on 
the incompleteness related to quantum gravity and show 
that this limitation can be overcome. In addition, we 
find that pre-inflationary dynamics can produce certain 
deviations from the BD vacuum at the onset of inflation, 
leading to novel effects which could be seen, e.g., in non- 
Gaussianities through future measurements of the halo 
bias and the '/x-type distortions' in the CMB [2]. 

Loop quantum gravity (LQG) offers a natural frame- 
work to address these issues because effects of its under- 
lying quantum geometry dominate at the Planck scale, 
leading to singularity resolution in a variety of cosmolog- 
ical models, including some that admit anisotropies and 
inhomogeneities [3]. Even though LQG is still incom- 
plete, notable advances have occurred — e.g., in cosmol- 
ogy, analysis of black holes, and a derivation of the gravi- 
ton propagator — by using the following strategy: First 
carry out a truncation of the classical theory geared to the 
given physical problem and then use LQG techniques to 
construct the quantum theory [4]. For inflation, then, we 
are led to focus just on first order perturbations off the 
spatially flat Friedman backgrounds with a scalar field 
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<t>. In numerical simulations we will use the quadratic 
potential V = (l/2)m 2 ^ 2 with m = 1.21 x 10 _6 mpi, the 
value that comes from the 7 year WMAP data [5, 6]. 
Throughout we use natural Planck units. 

The truncated phase space: We have P-prun = r o x Ti 
where r o is the 4-dimensional phase space of homoge- 
neous fields, and Ti, of the first order, purely inhomo- 
geneous perturbations thereon. Thus riy un is a normal 
bundle over r o . The base space r o is conveniently coor- 
dinatized by the scale factor a, the inflaton <j> and their 
conjugate momenta. Dynamics on r o is generated by the 
single, homogeneous, Hamiltonian constraint, C Q = 0. 
On Ti the first order constraints can be solved and one 
can readily pass to the reduced phase space Ti which we 
coordinatize by two tensor modes, collectively denoted 
by 7c in what follows, and the Mukhanov variable Qr 
representing the scalar mode. This passage to the re- 
duced phase space is carried out entirely in the phase 
space framework without using any equations of motion. 

Finally, a conceptually important point is that dynam- 
ics on TTrun, or on TTrun = r o x Fi, is not generated by 
a constraint, or, indeed by any Hamiltonian. Rather, 
the dynamical vector field X a on fxi-un has the form 
X a = fl^dpCo + Cl" p dpC 2 where Q a and are the 
symplectic structures on r o and F±, and <C' 2 is the part 
of the second order Hamiltonian constraint in which only 
terms that are quadratic in the first order perturbations 
are kept. X a fails to be Hamiltonian on Txrun because 
C2 depends not only on perturbations but also on back- 
ground quantities. However, given a dynamical trajec- 
tory 7o(t) on T and a perturbation at a point on it, X a 
provides a canonical lift of j {t) to the total space TTrun, 
describing the evolution of that perturbation along Yo(i). 

Quantum Kinematics: Since TTrun = r o x Ti, the to- 
tal Hilbert space is given by H = 1i ® Hi- The Hilbert 
space % of background fields consists of wave functions 
i b (a,4>) and its structure is well understood from loop 
quantum cosmology (LQC) [3]. For perturbations, we in- 
troduce an infrared cutoff so that A cu toff > A , the size of 
the observable universe. Physically, this amounts to 'ab- 
sorbing modes with A > A cu toff in the background'. Then 
there is a natural Hilbert space Hi on which perturba- 
tions Qr and 7^ act. It admits an infinite dimensional 
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sub-space of 4th order adiabatic states which are invari- 
ant under spatial translations, often called 'vacua'. Hi 
is generated by excitations on any one of them. (For an 
alternate characterization see [7].) Note however that, in 
contrast to quantum field theory on strictly stationary 
space-times, Hi does not have a preferred vacuum state, 
or a canonical notion of particles. 

The key difference from standard inflation is that quan- 
tum fields Qk , Tk now propagate on a quantum geometry 
represented by * (a, <P) rather than on a classical Fried- 
mann solution (a(t), 4>(t)). These quantum geometries 
are all regular, free of singularities. Thus, by construc- 
tion, the framework encompasses the Planck regime. 

Now, the quantum geometry underlying LQG is sub- 
tle [4]. For example, while there is a minimum non-zero 
eigenvalue of the area operator, there is no such mini- 
mum for the volume operator, although its eigenvalues 
are also discrete. In the present truncated theory, per- 
turbative modes with arbitrarily high frequencies are al- 
lowed even though there is a quantum geometry f in 
the background. By itself, this is not a problem. In 
our homogeneous sector, for example, the inflaton mo- 
mentum can be arbitrarily large but still the energy 
density p is bounded above by p 

max l ~^ J 

0.41/9 P1 [3]. The 

real trans-Planckian issue for us is whether the energy 
density in perturbations remains (not only bounded but) 
small compared to the background all the way back to 
the bounce. Only then would we be assured of a self con- 
sistent solution, justifying our truncation which ignores 
the back reaction. Otherwise one would have to await a 
full quantum gravity theory. 

Quantum dynamics: Since the classical dynamics on 
TTrun is not generated by a constraint, contrary to what 
is often done, one cannot recover quantum dynamics for 
the total system by imposing a quantum constraint. As 
in the classical theory, we can do this only in the homo- 
geneous sector H Q and we then have to 'lift' the resulting 
quantum trajectory to the full %. On H. one can follow 
the standard procedure in LQC. It again leads us to rein- 
terpret the quantum Hamiltonian constraint C Q f = 
as an 'evolution' equation, —ihd < p 1 i' (a,(j)) = H a ^ (a,<j)), 
with respect to the relational or emergent time variable 
4> generated by a time dependent Hamiltonian H Q [3, 8]. 
However, relative to the simpler massless case [3], one en- 
counters certain technical complications because of the 
presence of a potential. Their origin and resolution is 
analogous to that in the case where <fi is massless but 
there is a positive cosmological constant [9] . 

In the V(4>) — case, detailed investigations have 
shown that wave functions ^ (a, <f>) of physical interest 
remain sharply peaked even in the Planck era and fol- 
low quantum corrected effective trajectories. For V(<f>) — 
(l/2)m 2 2 , solutions to effective equations continue to 
undergo a bounce when p = /9 max and to agree with gen- 
eral relativity for p < 10~ 3 p max . However, because of 
computational limitations, so far the quantum wave func- 
tions \l/ (a, </>) have been calculated only when the bounce 
is kinetic energy dominated [8]. In this case, the peaks 



of wave functions of interest again follow the effective 
trajectories as expected. We restrict ourselves to back- 
ground quantum geometries *f/ (a,(j>) with this property. 
Each of them provides a probability amplitude for var- 
ious classical space-time geometries to occur. They are 
peaked not on classical Friedmann solutions but rather 
on quantum corrected bouncing solutions. Furthermore, 
there are fluctuations around these peaks. The challenge 
is to capture the effects of this background quantum ge- 
ometry ^> (a,<p) on the dynamics of perturbations. 

To meet it, we use the conceptual framework of quan- 
tum field theory on cosmological quantum geometries, 
introduced in [10]. An extension of that framework to in- 
corporate an infinite number of modes, with appropriate 
regularization and renormalization, provides the dynam- 
ical equation for states ip(Q%, T£) of perturbations on the 
background quantum geometry ^> . A key result is that 
this evolution is equivalent to that of test perturbations 
propagating on a dressed, effective, smooth metric 

g ab dx a dx b = ds 2 = d 2 ((f)) (-dfj 2 + dx 2 ) 

where the dressed scale factor a and the dressed confor- 
mal time fj are given by 

a 4 = (H~ h a 4 (0) HohiHoY 1 ; dfj = d 2 (<f>) (H- 1 ) dcf>. 

This result is exact within our truncation scheme. It 
shows that the propagation of perturbations is sensitive 
to properties of the state <I/ even beyond the quantum 
corrected effective geometry followed by its peak; it also 
senses quantum fluctuations around this peak. However, 
interestingly, this dependence is neatly coded in just two 
'dressed' quantities, rj and a. This is analogous to the fact 
that although light propagating in a medium interacts 
with its atoms, the net effect can be captured in just a 
few parameters such as the refractive index. 

This result greatly simplifies our task conceptually and 
enables us to use the technical tools of mode by mode 
regularization and renormalization from the well devel- 
oped adiabatic scheme of quantum field theory on classi- 
cal cosmological space-times [11]. For tensor modes, for 
example, one obtains the following evolution equation 



^fd 3 k %\h\ 2 + 



i°L \fA 2 

4k I'fcl 



C k (fi) MTfrij) 



where p£ is the momentum conjugate to 7%, n — 8nG 
and Ck{fj) arc c- numbers, derived from the 4th order adi- 
abatic regularization that depend only on k = \k\. 

Initial conditions: Since the big bang is replaced by 
the big bounce, it is natural to specify initial conditions 
at the bounce. The initial state can be taken to be of 
the form \I/ ® ip because perturbations are treated as 
test fields. This tensor product form prevails so long as 
the the back reaction remains negligible during evolution. 
To specify the initial condition for ^> a let us first recall 
that, in effective LQC, all dynamical trajectories enter a 
slow roll phase compatible with the 7 year WMAP data 
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FIG. 1: Ratio of our LQG power spectrum for scalar per- 
turbations to the standard inflationary power spectrum. For 
small k, the ratio oscillates rapidly with k. The solid curve 
shows averages over bins of width Afc = 0.5 ^pi" 1 . The inset 
shows a blow-up of the interesting region around k = 9. 



unless 0b , the value of the infiaton at the bounce, lies in a 
very small region R of the constraint surface [6] . We will 
assume that, at the bounce, the background quantum 
state $ o is sharply peaked at a point on the constraint 
surface outside this R. In this sense the initial data for 
\l/ is generic. For perturbations, we assume that the 
initial tp is a 4th order adiabatic 'vacuum' such that the 
expectation value of the renormalized energy density in 
tp is negligible compared to that in the background. This 
is a large class of initial data for test fields, selected by 
general symmetry requirements. 

Physically, we are assuming 'initial quantum homo- 
geneity' i.e., requiring that the region which expands to 
become the observable universe is homogeneous at the 
bounce except for 'vacuum fluctuations'. While this is a 
strong restriction, it may be naturally realized in LQG 
because: i) In solutions of interest, the observable uni- 
verse has a radius < 10£p\ at the bounce; and, ii) The 
strong repulsive force due to quantum geometry that 
causes the bounce has a 'diluting effect'. It could make 
this 'quantum homogeneity' generic, 'washing out' the 
memory of the pre-bounce dynamics at the scale < 10£pi. 

Our remaining task is two folds: i) starting from these 
initial conditions, calculate the power spectrum for scalar 
and tensor modes at the end of the slow roll inflation; 
and, ii) verify if the back reaction continues to remain 
negligible all the way to the onset of the slow roll so that 
our initial truncation is a self consistent approximation. 

Power spectrum: As noted above, in bounces with ki- 
netic energy domination on which we focus, the quan- 
tum state ty (a,4>) is known to remain sharply peaked 
on effective trajectories. Therefore, in numerical simula- 
tions a 'mean field' approximation was made by replacing 
5(0), r){4>) by the mean values of these operators. For the 
background, several simulations were carried out with 4>b 
in (0.93mpi, 1.5mpi), which, as we will see below, is the 
most interesting range. Three sets of initial conditions 
were chosen for the quantum state tp of perturbations. 



The power spectrum at the end of inflation was computed 
in each case for both scalar and tensor modes. Results 
are all very similar. FIG.l shows how the LQC scalar 
power spectrum relates to the prediction of standard in- 
flation for the case where 4>b — 1.15mpi, and the initial 
state tp is the 'obvious' or 'standard' 4th order adiabatic 
vacuum. We found that the plot is largely insensitive to 
choices of initial conditions within our class. 

Recall, however, that the 7 year WMAP data [5] covers 
only a window (fc m ; n w £;*/8.58, A; max « 2000fc m i n ) in the 
co-moving k space. Here the reference mode k* is the one 
that exits the Hubble radius at time fjk* when the Hub- 
ble parameter is given by H(rjk*) — 7.83 x 10~ 6 mpi. In 
FIG.l, numerical values of the co-moving k were calcu- 
lated using the scale factor convention ob = 1, rather than 
atoday = 1- (The physical wave numbers are of course 
convention independent). In each simulation, we first 
locate the scale factor a(fjk*) by setting H = H(fjk*), 
and then determine k* via k* — a(rjk*)H(fj) c *). Since we 
have ob — 1, values of a(fjk*) and k* depend on the pre- 
inflationary background dynamics which turns out to be 
governed entirely by 4>B- Therefore, in FIG.l the obser- 
vationally relevant window depends on the value of (pB , 
moving steadily to right as 0b increases. 

The plot has two interesting features. First, the LQG 
power spectrum is virtually indistinguishable from that of 
standard inflation if k m - m > 9mpi. This occurs when 
<Pb > 1.2mpi. Second, for smaller values of fc m in, the 
observational window admits modes for which the two 
power spectra are noticeably different. For concreteness, 
let us set cpB — 1.15mpi. Then A: m i n ~ 1.07mpi and these 
modes correspond to £ < 30 in the WMAP angular de- 
composition for which observational error bars are large. 
Therefore the LQG power spectrum is also viable but the 
predicted quantum state of perturbations at the onset of 
inflation is not the BD vacuum for <pB < 1.2mpi. 

Self consistency: Whether the test field approximation 
continues to hold in the Planck regime is an intricate 
issue and had not been explored before. FIG. 2 shows 
that there do exist solutions ip in which the renormal- 
ized energy density in perturbations remains low com- 
pared to the background throughout the evolution. (Here, 
we have set fc C utoff = k = 30mpi, which corresponds to 
(pB ~ 1.23topi.) Furthermore, there is an analytical ar- 
gument showing that every such tp admits a well-defined 
neighborhood (in the infinite dimensional space of the 4th 
order adiabatic 'vacua') with the same property. This es- 
tablishes that if cpB > 1.23mpi, the truncated theory ad- 
mits a rich set of self consistent solutions. Furthermore, 
in each of them tp has extremely small excitations over 
the BD vacuum at the onset of slow roll. These solutions 
provide viable extensions of the standard inflationary sce- 
nario all the way to the Planck scale. 

What about the small but interesting window (pB < 
1.2mpi? So far we only have upper bounds for the renor- 
malized energy density in perturbations and these are far 
from being optimal. Therefore we do not yet have an ex- 
plicit solution establishing the validity of the test field 
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FIG. 2: Energy density in the background (upper curve) and 
an upper bound on the renormalized energy density in per- 
turbations (lower curve) are plotted against time from the 
bounce to the onset of slow roll, using Planck units. The test 
field approximation holds across a change of over 11 orders of 
magnitude in both quantities. 



approximation in this window. 

Summary and Discussion: Using LQG ideas and tech- 
niques, we have extended the inflationary paradigm all 
the way to the deep Planck regime. At the big bounce, 
one can specify natural initial conditions for the quan- 
tum state ^f a that encodes the background homogeneous 
quantum geometry, as well as for ip that describes the 
quantum state of perturbations. There is a precise sense 
in which generic initial conditions for the background 
lead to a slow roll phase compatible with the 7 year 
WMAP data [6] . We have shown that there is a large set 
of initial data for -0 such that: i) at the onset of slow roll, 
ip is extremely close to the BD vacuum, and, ii) the test 
field approximation behind the truncation strategy holds. 
Each of these solutions provides a viable quantum gravity 
completion of the standard inflationary paradigm. How- 
ever particle physics issues still remain. 

In addition, there exists a narrow window in </>b for 
which the quantum state ip at the onset of inflation has 
an appreciable number of 'BD particles' (but within the 
current observational limits) . The physical origin of this 



effect can be explained in terms of the new scale fcjvf de- 
fined by the universal value of the scalar curvature at 
the bounce. Excitations with k < fcjvf are created in 
the Planck regime near the bounce. It turns out that if 
the number N of e-foldings in a between the bounce and 
V = Vk* is less than 15, then fc m i n < whence some of 
these modes would be in the window accessible through 
CMB. N < 15 corresponds to <pB < 1.2mpi, precisely 
the regime in which the LQC power spectrum is different 
from the BD vacuum. Future measurements should be 
sensitive to such deviations [2]. If they are observed at 
the scale k = fcjvf, the parameter space of initial condi- 
tions for ^ would be tightly squeezed, making much 
more detailed predictions feasible. In this sense, the 
framework expands the reach of observational cosmol- 
ogy all the way to the deep Planck regime. This general 
argument also shows that the pre-inflationary dynamics 
has negligible effect for modes with k ^> ku because their 
physical wave length turn out to be smaller than the cur- 
vature scale throughout the evolution. This explains the 
very close agreement between the LQC and the standard 
power spectrum at high k. 

Finally, interesting and complementary investigations 
of LQG dynamics between the bounce and the onset of 
slow roll have appeared in the literature recently (see, es- 
pecially, [12]). The distinguishing features of our analysis 
are: i) It is based on the general truncation strategy that 
has proven to be successful in other problems; ii) It pro- 
vides a systematic approach to quantum dynamics, made 
necessary by the fact that the classical evolution is not 
generated by a constraint on TTrun; iii) The treatment 
of initial states has been stream-lined; and, most impor- 
tantly, iv) While issues of regularization of the Hamil- 
tonian operator H\ , adiabatic renormalization of energy 
density, and consistency of the test field approximation 
were ignored so far, they have now been addressed using 
quantum field theory on quantum geometries. Details 
and subtleties which could not be included here will be 
discussed in two forthcoming articles. 
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